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Abstract. We have studied the supernova remnant G67.7+1.8 with the Chandra X-ray observatory. The remnant's X-ray 
morphology correlates well with the double-arc structure seen at radio wavelength. The X-ray spectra of the northern and 
southern rim of G67.7+1.8 exhibit emission line features of highly ionized metals, which suggests that most of the observed X- 
rays originate in a thermal plasma. We find magnesium , silicon, and sulphur overabundant relative to the solar values. Gaussian 
emission lines at ~ 4 keV and ~ 7 keV are detected. The ~ 4 keV line is consistent with K-emission lines from 44 Ca and/or 44 Sc 
whereas the ~ 7 keV line feature may arise from unresolved Fe-K lines. Chandra's sub-arcsecond angular resolution allowed 
us to detect four faint point sources located within ~ 1.5 arc-minutes of the geometrical remnant center. Among these objects, 
CXOU195424.75+312824.9 and CXOU195429.82+3 12834.1 do not have optical counterparts, leaving them as candidates for 
a possible compact stellar remnant. 
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1. Introduction 

Thanks to the progress made in high-energy astrophysics in 
recent years other manifestations of neutron stars besides 
accretion- and rotation-powered pulsars have been found. 
Among the various different categories known today, central 
compact objects (CCOs) are the most enigmatic ones. These 
sources show up as relatively faint X-ray point sources located 
close to the expansion center of their host supernova remnants 
(SNRs), strongly suggesting that they are indeed the compact 
remnants formed in the supernovae. Despite their apparently 
young age, the emission properties of CCOs are found to be 
very different from those observed in young rotation-powered 
pulsars. In the latter the X-ray emission is generally domi- 
nated by magnetospheric emission whereas the X-ray spec- 
tra of CCOs can be described well by composite blackbody 
models of temperatures (Ti,T<i) ~ (3 - 7) x 10 6 K and radii 
(R\,R.2) ~ 0.3 - 3 km of the projected blackbody emitting ar- 
eas that are much smaller than the canonical neutron star radius 
(see Becker 2008 for a review). In contrast to the young Crab- 
and Vela-like pulsars, no plerionic emission has been found 
around CCOs (see Hui & Becker 2006 and references therein). 

The temporal emission behavior of the known CCOs is an- 
other puzzle. So far, X-ray pulses have been found from three 
of the seven known sources. The detected periodic signals span 
a wide range from P ~ 0.1 s for the compact object in SNR Kes 
79 (Halpern et al. 2007) to P ~ 6.27 hrs in the case of the cen- 
tral point source in RCW 103 (De Luca et al. 2006). Different 
solutions have been proposed to explain their rotational dynam- 
ics. There is, however, no consensus so far on whether these 



objects e.g. were born with a slow spin and a weak magnetic 
field (Halpern et al. 2007) or at the other extreme that resem- 
bles the magnetars with a very short initial spin period and a 
very strong magnetic field (Li 2007). 

Currently, only seven SNRs have been confirmed as host- 
ing a CCO (cf. Tables 6.3 and 6.4 in Becker 2008). For a better 
understanding of their nature, the sample of CCOs needs to be 
increased. We therefore inspected known supernova remnants 
for centrally peaked X-ray emission by reanalyzing the ROS AT 
all-sky survey (RASS) data. G67.7+1.8 showed up in RASS 
data to have centrally bright X-ray emission that could poten- 
tially be associated with a compact source. A faint radio source 
seen in NRAO/VLA Sky Survey (NVSS) data (Condon et al. 
1998) close to the geometrical center of the remnant correlates 
well with the central X-ray emission (cf. Fig.l). The latter find- 
ing raised the question of whether this emission may eventually 
come from a young powerful pulsar and/or its synchrotron neb- 
ula. However, searching for a radio pulsar in G67.7+1.8 did 
not yield any detection. This non-detection placed an upper 
limit of 800 /vJy on pulsed radio emission at 600 MHz within 
G67.7+1.8 (Lorimer et al. 1998). Although this limit is not too 
constraining, it may suggest that the central emission could 
probably be associated with a CCO rather than with a young 
powerful rotation-powered pulsar. 

G67.7+1.8 was first discovered in a Galactic plane sur- 
vey at 327 MHz by Taylor et al. (1992), using the Westerbork 
Synthesis Radio Telescope. The radio image has shown that 
G67.7+1.8 has a double arc structure. The angular diameter is 
~ 9 . Its distance is quite uncertain. The only distance estimate 
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19 h 55 m s 54 m 30 s 54 m 3 53 m 30 s 
Right Ascension (J2000) 

Fig. 1. ROSAT PSPC RASS image of G67.7+1.8 (1. 0-2.0 keV) overlaid with contours from the NVSS survey data. The remnant 
appears center-filled in this data and has a bilateral radio arc-like structure. The peak of the central X-ray emission is found to be 
roughly matching the geometric center of the remnant. Faint radio sources are seen in the NVSS data near to the remnant center. 

available comes from the S-D relation, which puts the rem- the remnants distance might actually be smaller. In the same 
nant at ~ 17 kpc (Case & Bhattacharya 1998). However, optical work Mavromatakis et al. (2001) concludes that the energy re- 
filaments associated with the northern radio arc of G67.7+1.8 leased during the supernova explosion was significantly less 
were found by Mavromatakis et al. (2001) and indicated that than the canonical value of 10 51 ergs. 
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X-ray emission in the direction towards G67.7+1.8 was 
first observed in the ROSAT all-sky survey (Mavromatakis et 
al. 2001). Fitting the data with a thermal bremsstrahlung model 
indicated a plasma temperature in the range of ~ 0.2 - 0.3 keV. 

In this paper we present a detailed X-ray study of 
G67.7+1.8 and investigate the nature of its centrally peaked 
emission using Chandra ACIS-I observations. Its superior an- 
gular resolution enables us to perform a spectro-imaging study 
of the supernova remnant and to search for compact sources 
and/or a pulsar-wind nebula in its central part. 

2. Observations and data analysis 

G67.7+1.8 was observed with Chandra in 2007 March 4-5 
with the Advanced CCD Imaging Spectrometer (ACIS) using 
a frame time of 3.2 s. We utilized the whole ACIS-I CCD ar- 
ray to obtain an image of the remnant. The data were repro- 
cessed on 2007 August 7 in order to correct for a systematic 
aspect offset of ~ 0.4", possibly caused by a changing thermal 
environment. The total observing time on source was ~ 29 
ks. Examining the data for times of high background we no- 
ticed that the observation was affected by soft-proton flares. 
Cleaning the data by removing these flares from the data led to 
an effective exposure time of ~ 21 ks. The analysis was per- 
formed in the 0.5 - 8 keV energy band. 

2.1. Spatial analysis 

A false color image of the ACIS-I data of G67.7+1.8 is dis- 
played in Figure 2. Comparing the X-ray image with the radio 
contours obtained from the 1.4 GHz NVSS data (Condon et 
al. 1998) shows that the X-ray emission is entirely enclosed 
in the radio arc-like structure (see Figure 3). The center-filled 
X-ray emission and the radio shell-like morphology suggest 
G67.7+1.8 belongs to the category of "mixed-morphology" 
SNRs (Rho & Petre 1998). 

We have searched for the point sources in the whole ACIS- 
I data by means of the wavelet source detection algorithm 
and found 36 sources. These sources are marked by circles 
in Figure 2. The limiting count rate of the search was 2 x 
10~ 4 cts/s. The source positions, positional errors, signal-to- 
noise ratios as well as the estimates of extent of all these 36 
sources are given in Table 1. Cross-correlating these sources 
with the SIMBAD and NED databases did not result in any 
identification within a search radius of 5 arc sec around each 
source. 

Twelve of these newly detected sources are located within 
the supernova remnant. Four point-like sources, labeled in 
Figure 2 with the numbers #9, #10, #11 and #35 are lo- 
cated within ~ 1 .5 arcmin of the geometrical remnant center. 
Source #1 1 is the brightest among them. Hereafter, we desig- 
nate this source as CXOU195422.97+3 12902.1. It is also in- 
teresting to notice a trail-like radio feature close to the remnant 
center which has one end apparently coincide with source #10 
(cf. Figure 3). We designate it as CXOU195424.75+3 12824.9. 



1 See http://cxc.harvard.edu/cal/ASPECT/celmon/ for details 



Table 1. Positions of the X-ray sources in the field of ACIS-I 
as labeled in Figure 2. The corresponding Icr positional errors, 
signal-to-noise ratios as well as the estimates of source extent 
are also tabulated. 



Source 


RA (J2000) 


Dec (J2000) 


5 RA 


6 dec 


S/N" 


PSFRATIO " 




h:m:s 


d:m:s 


arc sec 


arcsec 






1 


19:54:57.172 


+31:28:36.59 


0.30 


0.22 


13.02 


1.08 


2 


19:54:53.728 


+31:32:12.81 


0.43 


0.20 


14.26 


1.47 


3 


19:54:52.507 


+31:24:33.74 


0.70 


0.28 


4.71 


0.80 


4 


19:54:46.466 


+31:31:26.13 


0.36 


0.47 


6.88 


2.38 


5 


19:54:40.352 


+31:23:11.84 


0.25 


0.26 


10.26 


0.75 


6 


19:54:37.457 


+31:32:56.14 


0.54 


0.33 


4.71 


2.19 


7 


19:54:35.356 


+31:23:49.09 


0.29 


0.19 


9.74 


0.83 


8 


19:54:28.210 


+31:32:37.34 


0.31 


0.32 


5.82 


2.76 


9 


19:54:24.920 


+31:30:04.42 


0.24 


0.16 


10.87 


5.02 


10 


19:54:24.745 


+31:28:24.96 


0.67 


0.40 


3.56 


9.29 


11 


19:54:22.965 


+31:29:02.08 


0.14 


0.17 


21.18 


6.13 


12 


19:54:21.016 


+31:33:33.80 


0.13 


0.12 


24.81 


1.35 


13 


19:54:16.717 


+31:32:11.63 


0.33 


0.28 


6.12 


2.15 


14 


19:54:12.784 


+31:24:30.48 


0.37 


0.52 


8.36 


1.78 


15 


19:54:10.148 


+31:30:01.07 


0.51 


0.37 


2.67 


1.99 


16 


19:54:09.741 


+31:24:11.62 


0.20 


0.39 


3.45 


0.47 


17 


19:54:02.592 


+31:27:46.74 


0.27 


0.34 


4.59 


1.03 


18 


19:54:00.649 


+31:23:28.01 


0.26 


0.34 


16.53 


0.95 


19 


19:53:58.504 


+31:26:19.25 


0.35 


0.50 


9.68 


1.31 


20 


19:53:52.939 


+31:25:13.26 


0.28 


0.40 


15.18 


0.90 


21 


19:53:52.315 


+31:27:08.97 


0.34 


0.40 


13.01 


1.22 


22 


19:55:07.125 


+31:30:03.42 


0.85 


0.70 


2.89 


0.81 


23 


19:55:03.136 


+31:26:12.22 


0.82 


0.51 


3.28 


0.96 


24 


19:54:57.325 


+31:28:09.59 


0.55 


0.37 


6.97 


1.36 


25 


19:54:49.830 


+31:26:45.24 


0.61 


0.38 


3.19 


1.20 


26 


19:54:30.035 


+31:31:52.60 


0.73 


0.44 


4.22 


5.44 


27 


19:54:20.292 


+31:21:33.27 


0.53 


0.54 


4.96 


0.85 


28 


19:54:19.929 


+31:22:50.35 


0.66 


0.36 


4.88 


1.28 


29 


19:54:06.334 


+31:27:18.46 


0.42 


0.43 


4.03 


1.68 


30 


19:53:58.555 


+31:31:33.96 


0.45 


0.60 


4.12 


1.19 


31 


19:53:58.106 


+31:32:01.81 


0.53 


0.40 


4.25 


1 .09 


32 


19:53:53.202 


+31:24:43.93 


0.60 


0.57 


4.45 


0.68 


33 


19:53:47.918 


+31:28:22.49 


0.81 


0.60 


3.66 


1.12 


34 


19:55:03.855 


+31:29:04.65 


0.54 


0.63 


4.10 


0.90 


35 


19:54:29.821 


+31:28:34.14 


0.67 


0.49 


3.78 


10.81 


36 


19:54:15.732 


+31:33:27.84 


0.50 


0.42 


4.38 


1.56 



" Estimates of source significance in units of Gehrels error: 

cr G = 1 + VCg + 0.75 where Cg is the background counts. 

b The ratios between the source extents and the estimates of the PSF 

sizes. 



This source appears to be slightly extended, which is also sug- 
gested by the ratio between its extent and the estimated PSF 
size (cf. Table 1). However, the low significance level of the 
detection and the patchy envirnoment of the remnant emission 
do not aUow one to determine if CXOU195424.75+3 12824.9 
is intrinsically extended. In the false colour image, source #35 
appears to have the hardest X-ray emission among these four 
central point sources (see Figure 2). We designate this source 
as CXOU195429.82+312834.1 hereafter. 

2.2. Spectral analysis of G67.7+1.8 

Before we have extracted the remnant spectra, all point-like 
sources were removed. The spectra were then extracted from 
the elliptical-shaped regions illustrated in Figure 4. These re- 
gions cover both the northern and the southern radio shells. 
The background spectra for the corresponding CCD chips were 
extracted from the boxed regions marked in Figure 4. We uti- 
lized the tool SPECEXTRACT in the data reduction package 
CIAO 3.4.1 with the calibration data CALDB 3.4.1 to extract 



4 



C. Y. Hui and W. Becker: X-ray observations of G67.7+1.8 




Fig. 2. Chandra ACIS-I false color image of G67.7+1 .8 (red: 0.3-0.75 keV, green: 0.75-2 keV, blue: 2-8 keV). The binning factor 
of this image is 2". It was adaptively smoothed by a Gaussian kernel of cr < 6". 36 point sources are detected in the field-of-view 
(cf. Table 1). Source #1 1 is designated as CXOU195422.97 +312902.1, source #10 is CXOU195424.75+3 12824.9. Top is north 
and left is east. 



the spectra and to compute the response files. After background 
subtraction there are 1028 and 610 net counts available for the 
spectral analysis of the northern rim (i.e. the sum of two north- 
ern elliptical regions) and southern rim (i.e. the sum of two 
southern elliptical regions), respectively. The background con- 
tributions are found to be ~ 60% in both, the northern and 
southern rims. We binned all the remnant spectra dynamically 
so as to have at least 50 counts per bin. All the spectral fit- 
tings were performed with XSPEC 11.3.2. The parameters of 



the best-fit models are summarized in Table 2. All quoted errors 
are lcr for 1 parameter of interest. 

2.2.1 . The northern rim 

The X-ray spectrum of the northern rim is shown in Figure 5 
which appears to be rich with line features. We attempted to 
fit the data with an absorbed collisional ionization equilib- 
rium (CIE) plasma model (XSPEC model: VEQUIL). To ex- 
amine the abundance of metals, we thawed the parameters 
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Fig. 3. 15x15 arcmin field of G67.7+1.8 as observed by Chandra ACIS-I in 0.5-8 keV. The radio contours at the levels of 0.6-6 
mJy/beam are from the NVSS data. The blue circles indicate the locations of four point sources near the center of G67.7+1.8. 
The inset displays the raw image of 4 x 4 arcmin of the remnant center with the X-ray sources labeled accordingly. 



individually to see whether the goodness-of-fit can be im- 
proved and/or the abundance significantly deviates from the so- 
lar value. However, this single component model cannot depict 
the data beyond 3 keV (with ^ 2 =45.32 for 33 D.O.F.) and re- 
quires an unreasonably overabundance of iron and calcium (i.e. 
hundred times of the solar values). 

We proceeded to examine whether the excess in the residual 
can be modeled by a second thermal component. Adding an- 
other CIE component with the solar abundances, we found no 
improvement in the goodness-of-fit (^ 2 =47.71 for 31 D.O.F). 

After examining the residual carefully, we found that com- 
bining a CIE model with two additional Gaussian components 
can model the observed spectrum. The composite model can 
describe the data within 0.5 - 8 keV very well: ;f 2 =25.41 for 29 
D.O.F. (cf. the lower panel of Figure 5). All the best-fit param- 
eters are tabulated in Table 2. 

The best-fit model yields a hydrogen column density of 
n H = (4.1 + 0.9) x 10 21 cirT 2 . Based on the Ha/H/3 line 
ratio, Mavromatakis et al. (2001) inferred an optical extinc- 



tion of ~ 2. This value implies a hydrogen column density 
of ~ 4 x 10 21 crrT 2 towards G67.7+1.8 (Predehl & Schmitt 
1995), which is in good agreement with our best-fit value. For 
comparison, the total galactic neutral hydrogen column den- 
sity towards G67.7+1.8 is ~ 10 22 cirT 2 (Kalberla et al. 2005; 
Dickey & Lockman 1990). The plasma temperature is found to 
be T - 6.6^Q 7 x 10 6 K. For the metal abundance, our analy- 
sis suggests that magnesium, silicon and sulphur are overabun- 
dant with respect to the solar values (Mg:2.6^J^, Si:2.8+J'2, 
S:13.6^ g). The best-fit parameters imply the unabsorbed flux 
to be f x = 6.7 x 10~ 13 ergs cirT 2 s" 1 in 0.5 - 8 keV. 

For the two additional Gaussian line features, the centroids 
of the line energy were found to locate at E\ = 4.0 + 0.2 keV 
and E% = 7.3^g 2 keV. The FWHMs of the line profiles are cr\ = 
0.3+° \ keV and <x 2 = 0.9+ 2 ° keV respectively. The best-fit line 
fluxes of the features are fy mei = SJ+j^xlO -6 photons crrT 2 s _1 
and /ii ne2 = 2.9+[ 2 4 7 x 10~ 5 photons crrT 2 s respectively. The 
possible physical nature and the significance of the emission 
line features are discussed in §3. 
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Fig. 4. Illustration of the elliptical-shaped regions used to 
extract the spectra from the northern and southern rim of 
G67.7+1.8. The background spectra were extracted from the 
rectangular regions. 

We have further examined whether a non-equilibrium ion- 
ization collisional plasma model (XSPEC model: VNEI) can 
provide a better fit than a CIE model. With a composite model 
of VNEI model plus two Gaussian components, We do not 
found any improvement in the goodness-of-fit (^ 2 =25.41 for 
28 D.O.F). All the best-fit parameters are consistent with those 
inferred from the CIE model. Moreover, the inferred ionization 
timescale is ~ 4.5 x 10 13 s crrT 3 which implies the condition 
of CIE (Borkowski et al. 2001). This justifies the validity of 
adopting the CIE model in our analysis. 

The detection of the emission line features prompted us to 
investigate their spatial distribution in G67.7+ 1 .8 by producing 
an X-ray image in the energy band of 3 -8 keV. However, owing 
to the high background contribution and the relatively limited 
photon statistics, we cannot identify any significantly excess 
emission from the background in the image. 

2.2.2. The southern rim 

The X-ray spectrum of the southern rim is shown in Figure 6. 
The spectrum appeared to be similar to what we observed from 
the northern rim, so that we fitted the data again with a CIE 
model. We have also examined the metal abundances as de- 
scribed in the previous subsection. However, the single compo- 
nent CIE model is not able to describe the data beyond ~ 4 keV 
(with^ 2 =38.82 for 23 D.O.F.). It also requires an unreasonably 
high abundance of iron. 

Inspecting the fit residuals reveals the existence of a 
Gaussian component centered at ~ 7 keV. However, modeling 
the data with the line energy of the Gaussian profile as a free 
parameter did not yield a stable fit. In particular, the centroid 




channel energy (keV) 

Fig. 5. X-ray spectrum of the emission from the northern rim 
(upper panel) observed by Chandra's ACIS-I2 and 13 CCDs. 
The residuals are from fitting an absorbed collisional ioniza- 
tion equilibrium plasma model plus two additional Gaussian 
components (lower panel). 

of the line cannot be constrained. This can be ascribed to the 
relatively low photon statistics of the southern rim spectrum. 
We therefore proceeded to fit the spectrum with the centroid 
of the line fixed at 7 keV. We found that the CIE model com- 
plemented with a Gaussian line at 7 keV can describe the data 
fairly well: ;^ 2 =31.37 for 23 D.O.F.. Different from the case of 
the northern rim, we did not find any excess around ~ 4 keV in 
the residual. 

The best-fit model yields a hydrogen column density of 
nu = 5.9^J j x 10 21 cirT 2 and a plasma temperature of T = 
(7.0 ± 0.7) x 10 6 K. Within the lcr errors these values are con- 
sistent with those inferred from the spectral fit of the north- 
ern rim. The overabundance of magnesium, silicon and sul- 
phur with respect to the solar values are Mg: 1 .7*1'), Si:4.7^, 
S:16.2^'' 4 7 . The best-fit parameters imply the unabsorbed flux 
to be /pi asma = 4.2 x 10~ 13 ergs cirT 2 s -1 in 0.5 - 8 keV. 

For the Gaussian line feature at 7 keV, the FWHM of the 
profile is found to be 1.1 * keV. The line flux of this feature is 
/line = (1.4 ± 0.5) x 10~ 5 photons crrT 2 s _1 . 

Adopting the VNEI model does not provide a better de- 
scription of the data than the CIE model. The composite model 
of VNEI component plus the Gaussian fixed at 7 keV re- 
sults in a goodness-of-fit of ^=31.36 for 22 D.O.F. All the 
best-fit parameters are consistent with those inferred from the 
CIE model. The ionization timescale is found to be ~ 4.5 X 
10 13 s cirT 3 which is similar to that inferred from the northern 
rim. Therefore, the CIE condition is justified and we will not 
further consider the VNEI model. 

2.3. Spectral analysis of the central point sources 

CXOU 195422.97+3 12902.1 (source #11) is the brightest one 
among the four point sources near to the remnant center. We 
extracted its spectrum from a circular region with a radius of 
10 arcsec centered at the source position (cf. Table 1). The 
background spectrum was extracted from a nearby region from 
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channel energy (keV) 

Fig. 6. X-ray spectrum of the emission from the southern rim 
(upper panel) observed by Chandra's ACIS-IO and II CCDs. 
The residuals are from fitting an absorbed collisional ionization 
equilibrium plasma model plus an additional Gaussian compo- 
nent with the centroid fixed at 7 keV (lower panel). 

a circle of 10 arcsec radius centered at RA=19 h 54 m 23.234 s , 
Dec=31°29'22.75" (J2000). After background subtraction, 54 
net counts were available for the spectral analysis of the source. 
Since the location of CXOU195422.97+3 12902.1 is close to 
the edge of ACIS-I3 CCD, we computed the response files 
manually with the CIAO tools MKARF and MKRMF. The 
spectrum was binned dynamically so as to have at least 5 
counts per bin. In view of the small photon statistic of these 
sources, we adopted the C-statistic (Cash 1979) for all the fit- 
tings. To better constrain the parameters, we fixed the hydro- 
gen column density at the optical extinction inferred value (i.e. 
4 x 10 21 cm" 2 ). 

Fitting with a power-law model results in a reasonable fit 
that yields a photon index of F = 2.7 ± 0.4 and a normalization 
at 1 keV of x 10~ 5 photons keV _1 cirr 2 s _1 . The best-fit 

column density is consistent with the values inferred from the 
spectra of the supernova remnant emission within the Icr error 
bound. The unabsorbed flux deduced for the best-fit power-law 
model parameters is f x = 4.5 x 10" 14 ergs crrr 2 s _1 in 0.5 - 8 
keV. 

A blackbody model can describe the spectrum of this 
source equally well. The best-fit model implies a temperature 
of T = 4.6^1 x 10 6 K. The radius of the projected blackbody 
emitting area is in the range ~ 1 70 - 310m and ~ 4 1 - 750 m 
for the adopted distances of 7 kpc and 17 kpc, respectively (see 
§3 for the discussion of the remnant distance). The unabsorbed 
flux deduced for the best-fit blackbody model is f x = 2.6x 10 -14 
ergs cirr 2 s _1 in 0.5 - 8 keV. 

For the other three fainter central X-ray sources it is in- 
teresting to compare their brightness and hardness with those 
of CXOU195422.97+3 12902.1. To do so, we have prepared 
their spectra and the response files in the same way as we did 
for CXOU195422.97+312902.1. Fixing the column density at 
nu — 4 x 10 21 cirr 2 we obtained the photon indices by fit- 
ting a power-law model to their spectra. The fitted parameters 
are summarised in Table 3. The photon index provides a mea- 



sure of the hardness of these X-ray sources. Whereas sources 
#9 and CXOU 195424.75 +3 12824.9 (source #10) are as soft 
as CXOU195422.97+312902.1, CXOU195429.82+3 12834.1 
(source #35) appears to show harder X-ray emission. We have 
also computed the absorption-corrected fluxes from the in- 
ferred power-law parameters which are given in Table 3. 

2.4. Optical identifications of the central point sources 

Given the limited photon statistics of the central point sources, 
the spectral analysis is not very constraining. Even for the 
brightest object CXOU195422.97+312902.1, we cannot dis- 
tinguish its emission nature between the thermal and the non- 
thermal scenarios. To investigate if these sources are promising 
neutron star candidates we proceeded to search for their optical 
counterparts by utilizing the USNO-B1.0 catalogue (Monet et 
al. 2003). 

For the search regions, we combine the positional errors of 
each source with the pointing uncertainty of the spacecraft. The 
uncertainty can be estimated from the distribution of aspect off- 
set for a sample of point sources with accurately known celes- 
tial position^ 2 ]. There is 68% of 70 sources imaged on ACIS-I 
have offsets smaller than ~ 0.4 arcsec. We adopted this value as 
the astrometric uncertainty and added to the quoted positional 
errors (cf. Table 1) in quadrature for each coordinate. 

For source #9, we have identified an optical counterpart at 
0.39 arcsec away from its X-ray position. It has a magnitude of 
B = 15.93 and R = 14.69 which implies an X-ray-to-optical 
flux ratio to be fx/ / opt ~ 10~ 3 . This ratio suggests source #9 is 
likely to be a field star which typically has a ratio fx/fopt < 0.3 
(Maccacaro et al. 1988). 

For CXOU195422.97+3 12902.1, an optical counterpart 
with B = 15.92 and R = 14.43 is consistent with its X-ray 
position. For this source the X-ray-to-optical flux ratio is at the 
level of ~ 10~ 3 which is in agreement with what is expected for 
a field star. 

We do not find any optical counterpart for 
CXOU195424.75+3 12824.9 and CXOU195429.82+3 12834.1. 
Adopting the limiting magnitude in the USNO-B1.0 catalogue 
to be 21 (cf. Monet et al. 2003), the upper limit of fx/fopt is 
found to be > 0.1. This upper limit is not particular constrain- 
ing. A deep optical observation is required to differentiate the 
nature of these two sources from the star with fx/fopt < 0.3 
and active galactic nuclei with fx/fopt < 50 (Maccacaro et al. 
1988; Stocke et al. 1991). 

3. Summary and Discussion 

We have performed a detailed spectro-imaging X-ray study of 
the supernova remnant G67.7+ 1 .8 with Chandra. Various prop- 
erties of G67.7+1.8, including distance, age and explosion en- 
ergy, are still poorly constrained. The type of the supernova 
and the nature of the progenitor are also unknown. With the 
Chandra observation it became possible for the first time to 
shed some more detailed light on the X-ray emission nature 
of G67.7+1.8 and its central point sources. 

2 http ://cxc . harvard. edu/cal/ASPECT/ celmon/ 
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Table 2. Spectral parameters of the best-fit model inferred from 
the Chandra observed spectra of G67.7+1.8. 







The Northern Rim 


The Southern Rim 


n„ (10 21 cm- 2 ) 

T (10 6 K) 

Mg 

Si 

S 

Norm cle (10- 4 )° 


4.1 ±0.9 

6 6 +0 5 
u.u _ 07 

2-8-11 
13.6% 
1 4+0.6 


5.9 + ! 2 

-1.1 

7.0 ±0.7 
1 7 +u 

1 • ' -0.7 

4 7+2.6 

^■'-1.7 

16 2+ 11 - 7 

1 o +0 - 5 


Ei (keV) 
a- 1 (keV) 
Norm gaussl 


(ltr 6 ) 6 


4.0 ±0.2 

3 +0 ' 2 

o q + 1.4 
-1.3 




E 2 (keV) 
o- 2 (keV) 

N0rm gau SS 2 


(io- 5 )* 


7 3 +32 
9 +2 ° 
2 9+ 12 - 7 


7 (fixed) 
1 l +04 
1.4 ±0.5 


x 2 

D.O.F. 




25.41 
29 


31.37 
23 



" The normalization of CIE model is expressed as (l(j-' 4 /42r£> 2 ) J N e N H dV where £> is the 
source distance in cm and N e and Ns are the post-shock electron and hydrogen densities 
in cm" 3 . 

h The normalization of the Gaussian model is in unit of photon cnT 2 s -1 . 



Table 3. Spectral characteristics of the central X-ray sources 
within the supernova remnant G67.7+1.8. 



Source 


net rate 
(10- 3 cts/s) 


r 


(ergs cm 2 


s- 1 ) 


9 


1.17 ±0.24 


3.2^;° 


2.8 x 10" 


■14 


10 


0.42 ± 0.26 


9 7+2.4 
^•'-1.7 


9.8 x 10" 


-15 


11 


2.54 ± 0.39 


2.7 ± 0.4 


4.5 x 10" 


■14 


35 


0.75 ± 0.23 


1 3 +0 - 7 

'■•■'-0.5 


1.2 x 10" 


■14 



The spectra of the remnants northern and southern rims can 
be described with a single temperature CIE model plus addi- 
tional Gaussian components. Our analysis suggests a plasma 
temperature of ~ 7 x 10 6 K. Assuming that the ions and elec- 
trons are fully equilibrated in the shock, the temperature im- 
plies that the blast wave has a velocity of ~ 800 km/s which 
is at the same order of magnitude as inferred from the RASS 
data (Mavromatakis et al. 2001). Although the single temper- 
ature plasma model describes our observed spectra very well, 
most supernova remnants require multi-temperature models to 
depict their X-ray spectra (e.g. Gaensler et al. 2008). A com- 
posite model can describe the shocked ejecta and the swept- 
up material from the surrounding with the components of high 
temperature and low temperature, respectively. However, for a 
more detailed modeling of the remnant spectrum, data with im- 
proved photon statistics are badly needed. 



Utilizing the Z-D relation, Case & Bhattacharya (1998) es- 
timated the distance to G67.7+1.8 to be ~ 17 kpc. However, the 
large uncertainties on the proportionality factor and the expo- 
nent of the 2 - D relation make the distance poorly constrained 
and results in a range of 7 - 27 kpc. On the other hand, the ob- 
servation of the optical filaments suggests a distance less than 
17 kpc, otherwise the detection of the optical radiation would 
be very difficult to explain in view of the measured interstellar 
extinction (Mavromatakis et al. 2001). In the following calcu- 
lations we assume the distance to be in the range of 7 - 17 kpc. 

The emission measure inferred from the spectral fits of the 
remnant emission allows us to estimate the hydrogen density, 
Nh, in the shocked regions. Assuming the post-shock densi- 
ties of hydrogen and electrons are uniform in the extraction 
region, the normalization of the CIE model can be approxi- 
mated by /V cie - 10~ w NnN e V/4nD 2 , where D is the distance to 
G67.7+1.8 in cm and V is the volume of interest in units of cm 3 . 
Assuming a geometry of oblated spheroid, the volume of inter- 
est for the northern rim and the southern rim are 3.9 x 10 55 <i, 3 

kpc 

cm 3 and 2.6 x 10 55 c/^ pc cm 3 , respectively. Here d^ pc is the rem- 
nant distance in units of kpc. Assuming a fully ionized plasma, 
NaN e can be approximated by 0.7A^. The best-fit normaliza- 
tions imply a post-shock hydrogen density of Na ~ 0.06 - 0.1 
cm- 3 for the assumed range of distances. 

Assuming G67.7+1.8 is in a Sedov phase, the shock tem- 
perature can be estimated by: 

r s -8.1xl0 6 (^J% 4 § K (1) 

where ti, E$\ and iV H _, are the time after the explosion in units 
of 10 4 years, the released kinetic energy in units of 10 51 ergs 
and the post-shock hydrogen density of 0.1 cm~ 3 respectively 
(McCray 1987). Mavromatakis et al. (2001) has argued the ex- 
plosion energy should be significantly less than the canonical 
value of 10 51 ergs. Adopting an explosion energy of E = 10 50 
ergs and the plasma temperature inferred from our spectral 
analysis as an estimate of T s (i.e. ~ 7 x 10 6 K), the range 
of N H ~ 0.06 - 0.1 cm" 3 yields the age of G67.7+1.8 to 
be ~ 5000 - 6000 yrs. On the other hand, an age bracket of 
~ 1 1000 - 13000 yrs is suggested from E = 10 51 ergs. 

Taking the estimates of the age, explosion energy and the 
density of the interstellar medium, we calculate the theoretical 
size of the remnant by (Culhane 1977): 

i 

/?,-2L5(J^-)%Jpc. (2) 

From the X-ray image, we estimated the angular radius of 
the remnant to be 9 S 5 arcmin. Comparing R s and 6 S , we can 
estimate the distance by d = R S /9 S . For an explosion energy 
of E = 10 50 ergs, the distance is estimated to be in a range of 
~ 7 - 12 kpc for N H = 0.06 - 0.1 cm" 3 and t = 5000 - 13000 
yrs. On the other hand, if £ = 10 51 ergs is adopted, it implies 
a distance bracket of ~ 12-20 kpc for the same ranges of Nn 
and t in the previous calculation. 

As aforementioned, the optical extinction suggests that the 
distance is at the lower side. On the basis of our estimation, 
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for a preferable distance much smaller than 17 kpc, we favor 
a scenario that the explosion is less energetic than a canoni- 
cal supernova (i.e. E < 10 51 ergs), which is consistent with 
the conclusion drawn by Mavromatakis et al. (2001) from their 
optical study. 

Assuming that the remnant G67.7+1.8 is the result of a low 
energy supernova explosion has interesting implications for the 
properties of its progenitor. In the context of the current under- 
standing in stellar evolution, a supernova with kinetic energy of 
~ 10 50 ergs can be a result of two different evolutionary tracks. 

In the first scenario, a Oxygen-Neon-Magnesium (ONeMg) 
core can be formed before Neon and subsequent nuclear burn- 
ing take place. If the neutrino cooling is efficient enough, the 
core temperature will be reduced and prevents further nuclear 
fusion. When the ONeMg core reaches the Chandrasakar mass, 
the electron degenerate pressure is no longer able to support 
the core. Furthermore, electron capture by 24 Mg and 20 Ne will 
further reduce the electron degenerate pressure in the core and 
trigger the core to collapse (Miyaji et al. 1980; Gutierrez, Canal 
& Garca-Berro 2005). This is known as the electron-capture su- 
pernova. Most simulations have shown that this type of explo- 
sion has a rather low energy (see Eldridge, Mattila, & Smartt 
2007). The progenitor's mass of an electron-capture supernova 
is limited in a narrow range of ~ 6-8 M (Eldridge, Mattila, & 
Smartt 2007). For more massive stars, they will go through all 
stages to silicon burning. In less massive stars, ONeMg cores 
cannot be formed. 

There is another possibility to produce a low energy super- 
nova. The exact evolution of a collapsing star after neutrino 
re-heating depends on the rate of early infall of stellar material 
on the collapsed core and on the binding energy of the enve- 
lope. If both are large, which is the case in high-mass stars, 
the energy required to accelerate and heat up the ejecta is not 
available, preventing a successful explosion and resulting in a 
supernova which appears to be under-energetic (see Heger et 
al. 2003; Eldridge & Tout 2004). The progenitor of this evolu- 
tionary track is likely to be > 20M o . 

We have also discovered emission line features in the rem- 
nant spectrum. The feature at E — 4.0±0.2 keV is detected only 
in the northern rim. The width of this feature (<x = 0.3 ke V) is a 
little larger than the energy resolution of ACIS-I at 4 keV (i.e. 
cr ~ 0.1 keV). This might indicate that the feature can be a 
blend of several lines. Its energy centroid is close to a number 
of K emission lines from 44 Ca and/or 44 Sc (cf. Table 5 in Iyudin 
et al. 2005 and references therein). Therefore, we speculate that 
the feature at 4 keV can possibly comprises these lines. 

If the identification is correct, it suggests a possible pres- 
ence of ^Ti because both 44 Ca and 44 Sc are produced in the 
decay chain 44 Ti— » 44 Sc— > 44 Ca. The half-life of 44 Ti is ~ 60 yrs. 
If the line feature at 4 keV is indeed from the decay of 44 Ti, its 
short half-life implies that the remnant should be rather young. 
Also, the production of 44 Ti in the supernova is very sensitive 
to the explosion mechanism and the ejecta dynamics (see dis- 
cussion in Iyudin et al. 2005). Therefore, obtaining the yield 
and spatial distribution of 44 Ti can help to further constrain the 
physical properties of the remnant as well as the nature of the 
progenitor. 



On the other hand, the broadening of the feature can also re- 
sult from the motion of the ejecta. However, this would imply 
the ejecta to have a velocity of ~ 10000 km/s. This is not con- 
sistent with G67.7+1.8 being the result of a low energy super- 
nova explosion as inferred from our X-ray data analysis as well 
from optical observation. Furthermore, the limited energy reso- 
lution of imaging data precludes any accurate measurement of 
the line width. High resolution grating spectroscopy is required 
to disentangle the possible blend of lines as well as determine 
their widths. 

Apart from the feature at 4 keV, we have also found an 
emission line feature at ~ 7 keV in both, the northern and 
southern rims. The large errors of the energy centroid and the 
width do not allow a firm identification. A deep observation is 
required to better constrain these parameters. We note that the 
approximate energy of the feature is rather close to the Fe-K 
emission. This leads us to speculate that the feature can be a 
result from unresolved Fe K a and lines. This also requires 
data with high spectroscopic resolution for further investiga- 
tion. 

The Chandra observation of G67.7+1.8 also leads to the 
detection of X-ray point sources at the proximity of the 
remnant's geometrical center. However, the limited photon 
statistics do not allow a constraining spectral analysis of 
these sources. The non-detection of optical counterparts for 
CXOU 195424.75 +3 12824.9 and CXOU195429.82+312834.1 
leaves these two objects as possible candiates for a stellar com- 
pact remnant. A dedicated deep optical observation is desirable 
to better constraining the nature of the point sources. 
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